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We use direct numerical simulations to compute structure functions, scaling exponents, probability 
density functions and turbulent transport coefficients of passive scalars in turbulent rotating helical 
and non-helical flows. We show that helicity affects the inertial range scaling of the velocity and 
of the passive scalar when rotation is present, with a spectral law consistent with ~ for the 

passive scalar variance spectrum. This scaling law is consistent with the phenomenological argument 
presented in [T] for rotating non-helical flows, wich states that if energy follows a E{k) ~ k~^ law, 
then the passive scalar variance follows a law V{k) ~ with ng = (5 — n)/2. With the second 
order scaling exponent obtained from this law, and using the Kraichnan model, we obtain anomalous 
scaling exponents for the passive scalar that are in good agreement with the numerical results. 
Intermittency of the passive scalar is found to be stronger than in the non-helical rotating case, 
a result that is also conhrmed by stronger non-Gaussian tails in the probability density functions 
of field increments. Finally, Pick’s law is used to compute the effective diffusion coefficients in the 
directions parallel and perpendicular to the rotation axis. Calculations indicate that horizontal 
diffusion decreases in the presence of helicity in rotating flows, while vertical diffusion increases. 

We use a mean field argument to explain this behavior in terms of the amplitude of velocity field 
fluctuations. 


I. INTRODUCTION 

The study of passive scalar advection, mixing and dif¬ 
fusion by anisotropic turbulence has gained more and 
more relevance over the years. Nowadays, it is well 
known that passive scalars share similarities with three- 
dimensional Navier-Stokes turbulence lai, presenting 
a direct cascade, anomalous scaling and intermittency 
nil]. Moreover, the study of passive scalar mixing in 
turbulent anisotropic flows is of interest in a wide vari¬ 
ety of geophysical and astrophysical problems, such as 
the transport of chemical elements in rotating stars 0- 
[5], the geodynamo |H], vertical transport and diffusion in 
the oceans nnun], and the transport of polutants and 
aerosols in the atmosphere |12) . 

Turbulent transport of passive scalars in rotating flows 
was previously studied in ffluani], although it has re¬ 
ceived less attention than the transport of passive scalars 
in isotropic turbulence [51I5HT5]. Moreover, the effect of 
helicity in the passive scalar transport in rotating flows 
has been practically ignored so far. It is known that 
helicity plays a key role in many problems such as in 
the dynamo effect mui], and the effect of flow helicity 
in the transport of passive vectors has been the subject 
of study in astrophysics for many years [18]. Results 
in niiiDi for isotropic turbulence indicate that passive 
scalar transport is sensitive to whether the flow is helical 
or not. In laminar flows, and in particular in biological 
flows, it has been found that helicity enhaces transport 
and mixing m- 

As helicity affects the direct cascade of energy in rotat¬ 
ing flows HD , leading to a steeper energy spectrum, it is 


to be expected that the passive scalar cascade to smaller 
scales should also be affected by the presence of helicity 
(see, e.g., jH])- From this point two questions naturally 
arise, which this work tries to answer: Is intermittency 
and the anomalous scaling of the passive scalar changed 
by the presence of helicity? And how is the transport 
and mixing of the passive scalar affected? While the for¬ 
mer question can be answered by computing scaling ex¬ 
ponents for rotating flows with and without helicity, the 
latter requires quantification of the turbulent transport 
in directions parallel and perpendicular to the rotation 
axis. 

The aim of this paper is then to characterize the turbu¬ 
lent scaling, transport and diffusion of passive scalars in 
rotating helical flows. To this end, we use data from di¬ 
rect numerical simulations of the Navier-Stokes equations 
in a rotating frame plus the advection-diffusion equation 
for a passive scalar. We use a spatial resolution of 512^ 
grid points in a regular periodic grid. 

The analisys is divided in two parts. First, to study 
the effect of helicity in the turbulent scaling laws of the 
passive scalar, we calculate velocity and passive scalar 
spectra. We compute structure functions for the velocity 
and the scalar using an axisymmetric decomposition, and 
consider the corresponding scaling exponents to quan¬ 
tify intermittency in each field. We also calculate prob¬ 
ability density functions (PDFs) for velocity held and 
passive scalar increments. As for non-helical rotating 
turbulence (see [1D)> ^Fat the passive scalar is 

more anisotropic than the velocity held at small scales. 
However, unlike the non-helical rotating case, the pas¬ 
sive scalar variance follows a spectral law consistent with 
~ where /c_l denotes wave vectors perpendicular 
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to the rotation axis. This scaling is shallower than the 
one found in the non-helical rotating case [23] . and is 
correctly predicted by a simple phenomenological rela¬ 
tion for the energy and passive scalar variance spectral 
indices. The passive scalar in the presence of helicity 
also becomes more intermittent than in the non-helical 
rotating case. 

Secondly, to study passive scalar diffusion, we com¬ 
pute effective anisotropic transport coefficients using the 
method used first in [21] for stratified flows, and later 
in [14] for rotating non-helical flows. Effective transport 
coefficients are obtained by studying the diffusion of an 
initial concentration of the passive scalar, and calculated 
using Pick’s law by measuring the average concentration 
and average spatial flux of the scalar as a function of time. 
For isotropic flows, we confirm that helicity increases tur¬ 
bulent diffusion (when compared with non-helical flows), 
in good agreement with previous studies and theoreti¬ 
cal predictions [H] [20] . In the presence of rotation, the 
overall effect of rotation (irrespectively of the content of 
helicity of the flow) is to decrease horizontal diffusion, 
while vertical diffusion remains approximately the same 
as in the isotropic case. Helicity further decreases hor¬ 
izontal diffusion, but slightly increases vertical diffusion 
(compared with the non-helical rotating case). The de¬ 
crease in horizontal diffusion is explained using a simple 
model for turbulence diffusivity based on the amplitude 
of the small scale velocity fluctuations. 


II. SET UP AND SIMULATIONS 

A. Equations and numerical method 

Data analyzed in the following sections stems form di¬ 
rect numerical simulations of the incompressible Navier- 
Stokes equations for the velocity u in a rotating frame, 
and of the advection-diffusion equation for the passive 
scalar 0, given by 


i9tu -I- u • Vu = —217 XU — Vp -I- vV'^VL + f, (1) 
V • u = 0, (2) 

dte + \i-^e = kV‘^9 -1- (j). (3) 

Here p is the pressure divided by the mass density (taken 
to be uniform in all simulations), v is the kinematic vis¬ 
cosity, and K is the scalar diffusivity. Also, f is an exter¬ 
nal force that drives the turbulence, (j) is the source of the 
scalar field, and O = Hz is the rotation angular velocity. 

To solve Eqs. 0-§ we use a parallel pseudospectral 
code in a three dimensional domain of linear size 27r 
with periodic boundary conditions [25[ [26]. The pres¬ 
sure is obtained by taking the divergence of Eq. (11, 
using the incompressibility condition given by Eq. (21, 


and solving the resulting Poisson equation. The equa¬ 
tions are evolved in time using a second order Runge- 
Kutta method. The code uses the 2/3-rule for dealias¬ 
ing, and as a result the maximum resolved wavenumber 
is kmax = A^/3, where N is the number of grid points in 
each direction. All simulations are well resolved, in the 
sense that the Kolmogorov dissipation wavenumbers for 
the kinetic energy and passive scalar variance, respectiv- 
elly ki, and k ^,, are smaller than the maximum wavenum¬ 
ber kmax at all times. More details of the numerical 
procedure can be found in [T]. 


B. Dimensionless numbers and parameters 


We will characterize the simulations using as dimen¬ 
sionless numbers the Reynolds, Peclet, and Rossby num¬ 
bers, defined as usual respectively as 


r. UL 

Re =- 

1 / 


( 4 ) 


Pe = —Re, 

hi 


( 5 ) 


Ro = 


U 


( 6 ) 


where U is the r.m.s. velocity, and L is the forcing scale 
of the flow defined as jkp with kp the forcing wave 
number. In all simulations U is close to unity in the 
turbulent steady state, and the kinematic viscosity is ^' = 
6 X 10“"^. The molecular scalar diffusivity is set equal to 
the kinematic viscosity for all runs, resulting in Pe = Re. 


C. Initial conditions and external forcing 

We performed a set of non-helical simulations and a 
set of helical simulations with varying Rossby numbers 
(see Table [I]) . In all cases, we first conducted a simula¬ 
tion solving only Eqs. Q and <§ (i.e., the incompress¬ 
ible Navier-Stokes equations without a passive scalar), 
starting from the fluid at rest (u = 0), and applying a 
random isotropic external mechanical forcing f to reach 
a turbulent steady state. This turbulent steady state 
was integrated for at least 13 turnover times. The me¬ 
chanical forcing f used to sustain the turbulent velocity 
field was a superposition of Fourier modes with random 
phases, delta-correlated in time, with tuneable injection 
of helicity using the methods described in Ref. EH- 

The procedure described above resulted in several runs 
as listed in Table [J with runs named with the letter “A” 
corresponding to simulations for wich the forcing inyected 
zero mean helicity, and runs labeled as “B” corresponding 
to runs with maximal injection of helicity. The final state 
of the velocity field in the turbulent steady state of these 
runs was used as initial condition for multiple runs in 
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TABLE I: Parameters used for the simulations: kp is the 
forcing wave number, Q is the rotation rate, Ro is the Rossby 
number, v is the kinematic viscosity. Re is the Reynolds num¬ 
ber, and = (u • V X u) is the mean helicity. Note that runs 
labeled with “A” have helicity fluctuating around zero, while 
runs labeled with “B” have non-zero helicity. 


Run 

kp 


Ro 

V 

Re 

H 

A1 

2 

0 

00 

6 X 10"^ 

525 

0 

A2 

2 

8 

0.02 

6 X 10"^ 

525 

0 

A3 

2 

16 

0.01 

6 X lO""* 

525 

0 

B1 

2 

0 

00 

6 X lO"'* 

525 

« 2 

B2 

2 

8 

0.02 

6 X lO""* 

525 

« 2 

B3 

2 

16 

0.01 

6 X lO”'* 

525 

« 2 


which the external mechanical forcing f was maintained, 
but a passive scalar was injected either as an initial con¬ 
centration 9(t = 0, x), or randomly injected in time using 
the source (j). 

These two different ways to inject the passive scalar 
depended on the properties of the scalar that were stud¬ 
ied. To characterize scaling laws and intermittency of the 
passive scalar in rotating helical and non-helical flows, 
the source term (j) was used to reach a turbulent steady 
state in the variance of the scalar as well as in the ki¬ 
netic energy. To this end, the source 4> was chosen as 
a superposition of Fourier modes with random phases, 
delta-correlated in time, injected at the same wavenum¬ 
bers kp used in the mechanical forcing f. 

Instead, to study passive scalar turbulent diffusion, 
and to compute effective transport coefficients, we turned 
off the source term in Eq. ([^ (i.e., we set (/) = 0). We then 
imposed two different initial conditions for the passive 
scalar, and integrated the velocity field and the passive 
scalar from those conditions to characterize horizontal 
and vertical diffusion. In each case, we used as initial 
condition Gaussian profiles as follows: 

0(t = O,x,) (7) 

where i = 1 or 3 (i.e., the initial profile can be a function 
solely of Xi = x, or solely oi x^ = z), = it (the pro¬ 

file is centered in the middle of the box, with the box of 
length 27r), and a = 1. When cci = x is used, this allows 
us to study the diffusion of the initial profile in the di¬ 
rection perpendicular to rotation (or “horizontal”), while 
when X 3 = z is used, we study diffusion in the direction 
parallel to rotation (or “vertical”). For a few runs, we 
verified explicitly that the diffusion in the x and y di¬ 
rections was the same (to be expected as rotating flows 
tend to be axisymmetric). These runs with no forcing 
and with Gaussian initial profiles for the scalar will be 
labeled with a subindex indicating the dependence of the 
initial profile (e.g., runs labeled Al^ or Al^ indicate the 
run A1 was continued with an initial Gaussian profile for 
0 that depends respectively on x or on z). 


III. TURBULENT SCALING LAWS 

In this section we present numerical results for the en¬ 
ergy and passive scalar spectra, structure functions, and 
PDFs for helical and non-helical rotating flows. To get 
the results in this section, the simulations in Table[l]were 
continued forcing the velocity and the passive scalar, to 
reach a turbulent steady state in both quantities. We 
first present the methods used to analyze the data, then 
presents the results for the spectra and inertial range scal¬ 
ing laws, and finally we characterize intermittency using 
structure functions and PDFs. We also compare the data 
with predictions from a simple phenomenological model, 
and from Kraichnan model for the passive scalar. 

A. Methods 

In this first part of the paper, the analisys consist on 
the characterization of flow anisotropy, scaling laws, and 
intermittency. To this end we consider power spectra, 
structure functions, and PDFs of passive scalar and ve¬ 
locity field increments for all runs. 

As a result of the anisotropy introduced by rotation, 
we consider reduced perpendicular energy and passive 
scalar spectra, namely E(kp) and V{k±_). These reduced 
spectra are defined by summing the power of all (velocity 
or passive scalar) modes in Fourier space over cylindrical 
shells with radius k_\_, with their axis aligned with the 
direction of the rotation axis. 

To compute structure functions and PDFs, field incre¬ 
ments must be defined first. Given the preferred direction 
introduced by rotation, it is natural to consider an ax¬ 
isymmetric decomposition for the increments. In general, 
the longitudinal increments of the velocity and the incre¬ 
ments of the passive scalar fields are defined respectively 
as 

(5u(x,l) = [u(x-bl) - u(x)] • |||, ( 8 ) 

(50(x,l) =0(x + l)-0(x), (9) 

where the increment 1 can point in any direction. Struc¬ 
ture functions of order p are then defined as 

5,(i) = (iMx,i)n, (10) 

for the velocity field, and as 

rp(i) = (|j0(x,i)n, (II) 

for the passive scalar field. Here, brackets denote spatial 
average over all values of x. 

These structure functions depend on the direction of 
the increment (i.e., they do not assume any symmetry 
in the flow). In simulations without rotation, the field is 
isotropic and the S'0(3) decomposition is used to calcu¬ 
late the isotropic component of the structure functions 
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]. In the rotating case, due to the axisymmetry of 
the flow, we will consider only increments perpendicular 
to z (the rotation axis), and increments parallel to z. We 
denote the former increments using l±, the latter with Zn, 
and we follow the procedure explained in detail in [TJ [3T] 
to average over several l± directions. 

This procedure to average Eqs. (10) and (11) over sev¬ 
eral directions can be summarize as follows. Velocity 
and passive scalar structure functions are computed from 
Eqs. ^ and Q using 26 different directions for the in¬ 
crements 1, generated by integer multiples of the vec- 
tors(l,0,0), (1,1,0), (2,1,0), (3,1,0), (0,1,0), (-1,1,0), 
(-1,2,0), (-2,1,0), (-1,2,0), (-1,3,0), (-3,1,0), 
(—1,3,0) (all vectors are in units of grid points in the sim¬ 
ulations), the 13 vectors obtained by multiplying them 
by —1, and the two vectors (0,0, ±1) for the transla¬ 
tions in z. Once all structure functions were calculated, 
the perpendicular structure functions Sp{l±) and Tp{l±) 
are obtained by averaging over the 24 directions in the 
X — y plane, and the parallel structure functions 5'p(Z||) 
and Tp{l\\) can be computed directly using the generators 
in the z direction. 

For all runs, this procedure was applied to Ns snap¬ 
shots of the velocity and of the passive scalar fields, sep¬ 
arated by at least one turnover time each. For large 
enough Reynolds number, the structure functions are ex¬ 
pected to show inertial range scaling, i.e., we expect that 
for some range of scales Sp ~ and Tp ^ Z^, where 
^p and (^p are, respectively, the scaling exponents of or¬ 
der p of the velocity and scalar fields. Scaling exponents 
shown bellow are calculated for all the snapshots ana¬ 
lyzed in each simulation, and averaged over time. Errors 
are then defined as the mean square error; e.g., for the 
passive scalar exponents, the error is 




iV, 


\ 


Ns 


E(Cp.-Cp)^ 


2=1 


( 12 ) 


where is the slope obtained from a least square fit for 
the Z-th snapshot, and (p is the mean value averaged over 
all snapshots. The error in the least square calculation 
of the slope for each snapshot is much smaller than this 
mean square error and neglected in the propagation of 
errors. Extended self-similarity [32l [33] is not used to 
obtain the scaling exponents. 


B. Energy and passive scalar spectra 

In the presence of rotation and in the absence of helic- 
ity, the spectral behavior of the passive scalar is strongly 
anisotropic and quasi-two dimentional [T] . As previously 
shown in P, E{k_L) - for the velocity field and 
V{k±) ^ for the passive scalar. The presence of 

helicity in rotating flows affects the cascade of energy and 
of the passive scalar to smaller scales. Numerical simula¬ 
tions in [22] showed that, when helicity is present in ro¬ 
tating flows, the direct cascade of helicity dominates over 




FIG. 1: (a) Reduced perpendicular passive scalar variance 
spectrum for run B3 (rotating and helical, solid) and for run 
A3 (rotating and non-helical, dashed), (b) Reduced perpen¬ 
dicular energy spectrum for run B3 (solid) and for run A3 
(dashed). 


the direct cascade of energy in the inertial range. This 
is the result of the development of an inverse cascade of 
energy, which leaves less energy available for the system 
to transfer to small scales. Assuming the direct cascade 
of helicity is dominant, a spectrum E{k±)H{k±) ^ kj^ 
is obtained from dimensional arguments [25]. In other 
words, if the energy spectrum satisfies E{k±) ~ ZcJ", 
then the helicity must follow a spectrum H{k±) ~ 

As a result, the energy spectrum becomes steeper as the 
flow becomes more helical, with the limit of a spectral in¬ 
dex n = 2.5 for the energy in the case of a turbulent flow 
with maximum helicity (in practice, this limit cannot be 
attained since in a flow with maximum helicity the non¬ 
linear term becomes negligible, resulting in no net energy 
transfer). 

Figure P shows the energy and passive scalar reduced 
perpendicular power spectra for runs A3 and B3 (both 
runs with rotation, and respectively without and with 
net helicity). The kinetic energy spectrum is steeper in 
the presence of helicity, compatible with E{kj_) ~ ZcJ^'^ 
scaling, while the passive scalar is close to V{k±) ~ 
scaling. Although resultion is moderate in these Simula- 
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FIG. 2: Reduced perpendicular energy spectrum compen¬ 
sated by k~^'^ (solid), reduced perpendicular passive scalar 
spectrum compensated by (dashed), and reduced per¬ 

pendicular helicity spectrum compensated by ® (dash- 
dotted) in run B3 (with rotation and helicity). 


tions (see [T] for more detailed studies of spectral scal¬ 
ing), the scaling laws can be further confirmed in Fig. 
where compensated energy and passive scalar spectra for 
run B3 are shown. In Fig. [^we also present the helic¬ 
ity spectrum compensated by (to confirm that the 

helicity spectral index and energy spectral index add to 
4). Similar scaling laws were observed in the rest of the 
helical rotating runs listed in Table |l] 

Following the phenomenological argument presented in 
[T] , we can explain the effect of helicity in the scaling of 
the passive scalar spectrum. From Eq. ([^, the passive 
scalar flux a can be estimated as 


a ~ 



(13) 


If we assume that the passive scalar has a direct cas¬ 
cade with constant flux a in the inertial range, then the 
passive scalar power spectrum V{k±) ^ df^/k^ can be 
estimated, using Eq. (13), as 


Ul, 


(14) 


For an energy spectrum E(k±) ^ and therefore for 
a characteristic velocity at scale I satisfying ui^ ~ 
the passive scalar spectrum in Eq. (14) results 


Vik±)^al] 


ak I 


(15) 


Therefore, the spectral index for the passive scalar iner¬ 
tial range is 


ng 


5 — n 
2 


(16) 


This simple phenomenological argument was proposed in 
[T] to take into account the effect of rotation in the spec¬ 
trum of the passive scalar. Here we confirm that this 



FIG. 3: Averaged axisymmetric structure functions (only for 
li_ increments) up to seventh order in run B3 (rotating and 
helical) for (a) the passive scalar, and (b) the velocity field. 


argument remains valid in the presence of helicity in the 
rotating flow. Moreover, when rotation is zero, we re¬ 
cover ng = 5/3, in good agreement with the Kolmogorov 
scaling previously observed for passive scalars in isotropic 
turbulence (see, e.g., la)- 


C. structure functions and scaling exponents 

Structure functions and scaling exponents for the pas¬ 
sive scalar in non-helical rotating flows were studied in 
detail in [T]. As a result, here we focus on the simu¬ 
lations with helical forcing. Figure shows the axisy- 
metric and perpendicular (i.e., only for perpendicular in¬ 
crements l±) structure functions for the passive scalar 
and for the velocity field up to seventh order for run B3. 
Each curve corresponds to an average over Ng = 8 snap¬ 
shots of the turbulent steady state of the simulation. The 
structure functions show a range of scales with approx¬ 
imately power-law scaling at intermediate scales, while 
at the smallest scales approach the ^ P scaling expected 
for a smooth field in the dissipative range. 

Figure shows a detail for the same run of the pas¬ 
sive scalar and velocity field second-order perpendicu¬ 
lar structure functions, respectively T 2 { 1 ±) and 82(1 ±), 
as well as the structure functions for increments par¬ 
allel to the rotation axis, T 2 { 1 \\) and S' 2 (/||). Stronger 
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P 


FIG. 5: Scaling exponents (with error bars) as a function 
of the order p in simulations of helical rotating turbulence, 
for the velocity held (triangles for run B2 and stars for run 
B3, both with helicity and with decreasing Rossby number), 
and for the passive scalar (diamonds for run B2 and squares 
for run B3). The solid line corresponds to the linear scaling 
expected for the velocity held exponents in the absence of 
intermittency, while the dash-dotted line corresponds to non- 
intermittent scaling for the passive scalar exponents. The dot¬ 
ted and dashed lines correspond to Kraichnan’s model with 
^2 = 0.4 and respectively with d = 2 and with d = 3. 


FIG. 4: Axisymmetric second order structure functions for 
run B3 (helical with rotation) for (a) the passive scalar, and 
(b) the velocity held. In both panels solid lines correspond to 
the parallel structure functions, while dashed lines correspond 
to the perpendicular structure functions. Slopes indicated 
as references correspond to the time average of the scaling 
exponents, obtained from a best ht in the inertial range of 
the structure functions at diherent times. 


anisotropy is observed at small scales for the passive 
scalar than for the velocity field, manifested as a larger 
difference between T 2 { 1 \\) and T 2 { 1 ±) than between 5'2(^||) 
and 82(1 ±)- Also, an inertial range with power-law scal¬ 
ing can be identified at intermediate scales in 12 (^j_) 
and 82(1 ±)- The range of scales is consistent with the 
wavenumbers of the inertial range in the corresponding 
spectra. The slopes indicated as a reference in Fig. 
correspond to the time average of the second-order scal¬ 
ing exponents, obtained from a best fit in the inertial 
range of all structure functions at different times. The 
second-order scaling exponents (in the perpendicular di¬ 
rection) are (2 = 0.41 ± 0.01 for the passive scalar, and 
^2 = 1.22 ± 0.01 for the velocity field. These values are 
in good agreement with the spectra V{k±) ~ and 
E{k±) ~ wich from dimensional analysis lead to 

T 2 { 1 ±} ~ and S'2 (/_l) ~ l];^. From the curves in 
Fig-i scaling exponents can also be computed for lower 
and higher orders. Based on the amount of statistics 
available, velocity and passive scalar exponents in the di¬ 
rect cascade range were computed for all runs up to the 
seventh order. 

Figure shows the resulting velocity scaling exponents 


Cp, and passive scalar exponents ^p, for runs B2 and B3 
(both helical, with Ro = 0.02 and 0.01 respectively). Lin¬ 
ear (non-intermittent) scalings for (^p and for are shown 
as a reference, based on the values of the second-order ex¬ 
ponents C2 and ^2- In Fig. we also show the the predic¬ 
tion of the Kraichnan model [4] , which is a model for the 
advection and diffusion of a passive scalar in a random, 
delta-correlated in time velocity field in a space with di¬ 
mensionality d. The scaling exponents for the passive 
scalar in this model are 


Cp — X '\/2dC2P + (d — C2)^ + (d—C2) 


2 1 


(17) 


For the curves in Fig. these exponents were evaluated 
with the value of ((2 obtained from the simulations, and 
using either d = 2 or d = 3. 

Scaling exponents for the velocity field are similar in 
both runs. The second-order velocity field exponent is 
^2 = 1.22 ± 0.02 for run B2, and 6 = 1-23 ± 0.01 for run 
B3. The velocity field exponents display the well-known 
deviations from linear scaling associated with intermit¬ 
tency, more evident for the higher order exponents and 
in the simulation with larger Rossby number (i.e., smaller 
rotation rate). The deviation from strict scale invariance 
is often quantified in terms of the intermittency exponent 
/i = 2^3 — ^6, which for these runs is /i = 0.6 ± 0.2 for run 
B2, and p = 0.2 ± 0.1 for run B3. The decrease in the 
values of p suggest a reduction of intermittency with in¬ 
creasing rotation, as observed before in simulations and 
in experiments [HinillMOH]- 

The passive scalar exponents for these two runs also 
display similar values.The second-order scaling exponent 
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FIG. 6: Scaling exponents (with error bars) as a function 
of the order p in simulations of rotating turbulence with and 
without helicity, for the velocity (triangles for run A3 without 
helicity, and stars for run B3 with helicity), and for the passive 
scalar (diamonds for run A3, and squares for run B3). 


is C 2 = 0.41 ± 0.01 for both runs. Deviations from linear 
scaling are observed, and the intermittency exponents are 
/is = 0.36 ± 0.06 for run B2, and = 0.27± 0.04 for run 
B3. For runs B2 and B3, Kraichnan’s model adjusts the 
numerical data best with C 2 = 0.4 and d = 2. The value 
of d is compatible with quasi-bidimensionalization in the 
spatial distribution of of the passive scalar in the presence 
of rotation, as reported in the presence of rotation in [T] . 

Overall, the decrease in the values of /i observed for 
both the velocity field and the passive scalar indicate a 
reduction of intermittency with decreasing Rossby num¬ 
ber. However, this reduction is more pronounced for the 
velocity field than for the passive scalar. 

Finally, we present a comparison between the scaling 
exponents in rotating turbulence with and without helic¬ 
ity. Figure shows the velocity field and passive scalar 
exponents for runs A3 and B3 (respectively without and 
with helicity). Deviations from linear (non-intermittent) 
scaling are larger for the passive scalar in B3, indicating 
stronger intermittency in the presence of helicity. 


D. Probability density functions 

Intermittency and small scale anisotropy can be also 
studied considering the PDFs of the field increments. In 
this section we present PDFs of longitudinal increments 
of the cc-component of the velocity field, as well as incre¬ 
ments and spatial derivatives of the passive scalar concen¬ 
tration. Quantities shown are normalized by their vari¬ 
ance, and a Gaussian curve with unit variance is shown 
as a reference. 

Figure shows the PDFs of the velocity and of the 
passive scalar increments for four different values of the 
spatial increment {I = 1.6, 0.8, 0.4, 0.2, and 0.1) in run 
B3 (with helicity). All the increments were considered 
in the x-direction (perpendicular to the axis of rotation, 
and for the velocity the cc-component was used to build 
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FIG. 7: Probability density functions in run B3, for five 
different horizontal spatial increments I — 1.6 (solid), 0.8 
(dashed), 0.4 (dash-dotted), 0.2 (dash-triple-dotted), and 0.1 
(long dashes), and for (a) the passive scalar, and (b) the x- 
component of the velocity field. A Gaussian curve with unit 
variance is indicated by the dotted curve. As intervals are 
decreased, curves depart more and more from the Gaussian 
distribution developing stronger tails. 



longitudinal increments. As a reference, and to compare 
with the values of the increments considered, the forc¬ 
ing scale in this runs is « tt, and the dissipative scale is 
« 0.05. Therefore, increments I = 0.8 and 0.4 correspond 
to scales in the inertial range. The PDFs of velocity and 
passive scalar increments for I = 1.6 are close to Gaus¬ 
sian, while for smaller spatial increments non-Gaussian 
tails develop. Note also that in the PDFs of passive scalar 
increments, a strong asymmetry develops for I = 0.4, 0.2 
and 0.1. 


IV. TURBULENT DIFFUSION 

In this second part of the paper, the aim is to char¬ 
acterize the turbulent diffusion of the passive scalar in 
rotating helical turbulence, and to compare it with turbu¬ 
lent diffusion in non-helical rotating turbulence, as well as 
with turbulent diffusion in isotropic turbulence. To this 
end, we simulate the flows starting from an initial Gaus¬ 
sian profile for the concentration of the passive scalar. 
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FIG. 8: (a) Averaged horizontal concentration 6 in run Ala; 
(no rotation, no helicity) at times t = Q (solid), 0.5 (dotted), 
1 (dashed), 1.25 (dash-dotted), and 1.5 (dash-triple-dotted), 
(b) Horizontal flux at the same times, (c) t) at the same 

times. 



FIG. 9: (a) Averaged horizontal concentration 6 in run Bla, 
(no rotation, helical) at times t = Q (solid), 0.5 (dotted), 
1 (dashed), 1.25 (dash-dotted), and 1.5 (dash-triple-dotted), 
(b) Horizontal flux at the same times, (c) 'D^ix, t) at the same 
times. 


and we let it diffuse in directions parallel and perpen¬ 
dicular to the rotation axis. We then quantify effective 
transport coefficients by measuring the time evolution of 
the averaged concentration, and using Pick’s law. 


A. Methods 

Before presenting the method used to measure the tur¬ 
bulent diffusion, we briefly recall how the simulations 
were conducted for this second study. As in the previous 
section, simulations in group A (see Table |T| correspond 


to simulations with zero mean helicity, while simulations 
in group B correspond to simulations with helical forcing 
and non-zero net helicity. As explained in Sec. |IV A[ for 
each run in the turbulent steady state of the velocity field, 
the simulation was extended twice with the same parame¬ 
ters and mechanical forcing, but with two different initial 
conditions for the passive scalar: a Gaussian profile for 
the concentration in the x-direction (to study horizon¬ 
tal diffusion), and a Gassial profile in the z-direction (to 
study vertical diffusion). To identifie these runs, an addi¬ 
tional subindex is used in this section to differentiate be¬ 
tween simulations with different dependence of the initial 
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FIG. 10: Horizontal turbulent diffusion as a function of time 
for runs Al^; (solid, no rotation and no helicity) and 
(dashed, no rotation but with helical forcing). 


Gaussian profile. As examples, a run labeled Al^, stands 
for a simulation with the parameters of run A1 (i.e., with 
zero mean helicity and no rotation) and with initial pro¬ 
file of the passive scalar in the cc-direction, while the label 
62^ indicates the run has helicity, rotation, and an initial 
dependence of the passive scalar in the z-direction. 

In each of these runs, we let the initial profile diffuse for 
several turnover times. Meanwhile, we compute and store 
quantities averaged over the two directions perpendicular 
to the direction over which the original Gaussian profile 
varies. In particular, we consider the averaged passive 
scalar concentration 6, and the spatial passive scalar flux 
0Ui, where z = I or 3 depending on the initial dependence 
of the Gaussian profile, and where the averages denoted 
by the overbars are done over the two remaining Carte¬ 
sian coordinates. Note the spatial flux Oui represents the 
amount of passive scalar transported in the z-direction 
per unit of time by the fluctuating (or turbulent) veloc¬ 
ity, since there is no mean flow in our simulations (we 
use delta-correlated in time random-forcing), Ui is the 
fluctuating velocity. 

Then, the pointwise effective turbulent diffusion coef¬ 
ficient is given by [2l] 


Vi{xi,t) = 


Ouj 

dx.O 


(18) 


This coefficient corresponds to how much passive scalar is 
transported by the fluctuating velocity per unit of varia¬ 
tion of 9 with respect to Xi. As already mentioned, z = 1 
stands for horizontal diffusion, while z = 3 stands for ver¬ 
tical diffusion, where the dependence on the direction of 
this coefficient is the result os the flow being anisotropic. 

From Pick’s law, the actual turbulent diffusion coeffi¬ 
cient is the average of 'Di{xi,t) over the coordinate Xi, 
and if the system is in a turbulent steady state, over 
time. From Eq. (181, we can define these averaged dif¬ 


fusion coefficients as follows. We can first average over 
the coordinate Xi to obtain a time dependent turbulent 




FIG. 11: (a) Averaged horizontal concentration 6 in run BSa,, 
at times t = 0 (solid), 0.25 (dotted), 0.5 (dashed), 0.75 (dash- 
dotted), and 1 (dash-triple-dotted), (b) Horizontal flux at the 
same times, (c) T>x(x,t) at the same times. 


diffusion. 


1 

^ y Vi{xi,t)Ax^, (19) 

and we can further average over time, to obtain the mean 
turbulent diffusion 

1 j-to+T 

P, = - / V,{t)dt. (20) 

J Jto 

Here, to and T are characteristic times of the flow. In 
practice, in our simulations the turbulent diffusion (t) 
first grows in time as the initial Gaussian prohle is mixed 
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FIG. 12: Horizontal turbulent diffusion as a function of time 
for runs A2x (solid) and B2a; (dashed) (Ro = 0.02, respec¬ 
tively without and with helicity). 

by the turbulence, then reaches an approximate steady 
state value for a few turnover times, and then decreases 
as the scalar becomes completely diluted (which happens 
after three or four turnover times). 

B. Isotropic helical turbulence 

In the absence of rotation, diffusion coeficients are ex¬ 
pected to be isotropic, and therefore horizontal and ver¬ 
tical turbulent diffusion should be the same within er¬ 
ror bars. Figure shows the mean passive scalar profile 
0{x,t), the horizontal flux 0ux{x,t), and the pointwise 
value of 'Dx{x,t), at five different times for run Al^, (no 
rotation and no net helicity). 

As time evolves, the mean profile 6{x,t) flattens and 
widens. The flux is antisymmetric: it is positive for x > ir 
and negative for x < tt. This behavior for the flux is to 
be expected, as at t = 0 there is an excess of passive 
scalar concentration at x = tt that must be transported 
by turbulent diffusion towards x = 0 and towards x = 2tt. 
The pointwise value of 'Dx{x, t) fluctuates around a mean 
value (which increases with time), except close to x = tt 
where it rapidly takes very large positive and negative 
values as in that point d^O approaches zero. The mean 
spatial value of V^ix^t) increases to its saturation value 
around to ~ 1.5; after this time it fluctuates around its 
value (see more details below). 

Figure shows the same quantities at five different 
times for run Bl^, (i.e., in a simulation without rotation 
but with injection of net helicity). The behavior of the 
mean concentration of the passive scalar, the horizontal 
scalar flux, and the pointwise value of V^^x, t) is qualita¬ 
tively the same as in the non-helical run Al^,. However, 
the helical run displays a larger diffusion of the mean con¬ 
centration of the scalar (as evidenced by the smaller max¬ 
imum value of 6[x,t) around x = tt and by the stronger 
tails close to to a; = 0 and 27r, when curves at the same 
time are compared in Figs. |^and[^. Also, the spatial 



FIG. 13: Horizontal turbulent diffusion as a function of time 
for runs ASa, (solid) and BSa, (dashed) (Ro = 0.01, respec¬ 
tively without and with helicity). 

flux Oux takes larger extreme values in the helical simu¬ 
lation, and the spatial average of V^ix, t) seems to result 
in larger values for the turbulent diffusion in this run. 

The increased turbulent diffusion in the presence of 
helicity is confirmed in Fig. |10[ which shows the hori¬ 
zontal turbulent diffusion as a function of time for runs 
Ala; and Bla;. In both runs T>x(t) grows from an ini¬ 
tially small value to its saturation value around to « 1.5. 
As observed above, turbulent diffusion saturates at sim¬ 
ilar times for the helical and the non-helical case, but to 
a larger value in the presence of helicity. Although in 
isotropic turbulence helicity does not affect significantly 
the energy scaling [22l ESI EH EQ] , an increase in the tur¬ 
bulent diffusion in the presence of helicity was predicted 
in |20j . Using renormalization group techniques, the au¬ 
thors estimated that turbulent diffusion in a helical flow 
can be up to a 50% larger than in a non-helical flow. In 
our simulations the averaged in time value of 'Dxit) is 
« 0.3 for run Alj,, and « 0.4 for run Bl^,, in reasonable 
agreement with the theoretical result. 

It is worth mentioning that the same analysis was per- 
fomed in simulations Al^ and BI 2 (i.e., the same runs 
but with an initial Gaussian profile in the z-direction). 
As expected from the flow isotropy, the same behavior 
was obtained. 


C. Rotating helical turbulence 

1. Horizontal diffusion 

Figure [TT] shows the mean profile of the passive scalar 
0{x,t), the horizontal flux 6ux{x,t), and the pointwise 
value of 'Dx{x,t) at five different times for run BS^, (heli¬ 
cal and with Ro = 0.01). In this case, note that the av¬ 
erage profile and the flux become asymmetric, i.e., there 
is an excess of concentration of 6{x,t) for x < tt, and 
the absolute value of the flux is larger tor x < tt than 
for X > TT. This assymetry is caused by the Coriolis force 
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FIG. 14: Averaged vertical concentration 6 in run BSz, at 
times t = 0 (solid), 0.25 (dotted), 0.5 (dashed), 0.75 (dash- 
dotted), and 1 (dash-triple-dotted), (b) Horizontal flux at the 
same times, (c) T>x{x,t) at the same times. 


and has been previously observed for rotating non-helical 
flows in [13 E]. In our runs, the passive scalar at t = 0 
is concentrated in a narrow band around x = it. The av¬ 
erage flux is thus towards positive values of x for x > tt 
and towards negative values of x for x < tt (i.e., in the 
direction of —V0, see, for instance, Fig. 16). The Coriolis 
force in Eq. 0 is — 2nz X u, and creates an overturning 
in in the x-y plane of the initially only dependent on x 
Gaussian profile, as will be shown later in more detail 
in spatial visualizations of the passive scalar. This over¬ 
turning also results in the asymmetry in Fig. 11 (for more 
details, see also [Ti]). 

By compiting the mean value of 'Dx{x, t) over the spa- 



t 

FIG. 15: Vertical turbulent diffusion as a function of time for 
runs A3z (solid) and B3z (dashed). The former run has no 
net helicity, while the latter has helical forcing. 


tial coordinate we obtain the turbulent diffusion coeffi¬ 
cient. Figure 12 shows first the horizontal turbulent dif¬ 
fusion as a function of time for runs A2x and B2a; (both 
with Ro = 0.02, without and with helicity respectively), 
and then Fig. shows the same quantity for runs AS^, 
and B3a; (Ro = 0.02, without and with helicity respec¬ 
tively). For both rotation rates, we observe that horizon¬ 
tal diffusion is smaller in the presence of helicity. This 
result is the opposite to that observed for the isotropic 
runs in the previous section, for which helicity increased 
the turbulent diffusion. 

As already mentioned, while in isotropic turbulence 
helicity does not affect the energy spectrum scaling (53 
123133 Eq] , in rotating turbulence the presence of helicity 
results in shallower horizontal spectrum for the energy, 
in comparison with rotating non-helical turbulence |221 
[23 • As a result, a smaller turbulent diffusion can be 
expected, as small scale velocity field flutuations should 
be less energeting in the helical rotating case. Indeed, in 
most two point closure models, the turbulent diffusivity is 
proportional to the mean kinetic energy in the turbulent 
fluctuations, u^/2, and if the kinetic energy spectrum is 
steeper, then the diffusivity should decrease. A simple 
mean field argument can illustrate this. We can split the 
velocity in a mean flow u, and a fluctuating component 
u', such that u = u -b u'. In our runs u = 0, and u = u'. 
Splitting the passive scalar in the same way we have 9 = 
9 + 0'. Replacing in Eq. |3l and averaging we obtain 


f)f) 

- = -V.(u^), (21) 

and subtracting this equation from Eq. 0 we then obtain 
r)f) 

-=-V.(u0). (22) 

We can integrate this last equation assuming the flow 
is correlated over the integral eddy turnover time r, to 
obtain 


9' « —rV • (u0) = —ru • Vd, 


(23) 
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FIG. 16: Passive scalar concentration in a horizontal slice of 
runs A3x (left) and B3x (right) at time t = 1. Note how the 
initial concentration (Gaussian, centered around x = n, and 
independent of the y-coordinate) gets distorted and diffussed. 


where incompressibility was used. 
Eq. 


Then, replacing in 



0 12 3 4 5 6 

y 



de d ,_, de 


(24) 


where the coefficient ruiUj can be interpreted as a tur¬ 
bulent diffusion. If the flow is isotropic, then V « ru^. A 
more refined mean field derivation of this expression can 
be found in [¥TH15] . while two point closure derivations 
can be found in [44l |43 . 

Although the argument above is only illustrative, it 
gives an interesting hint to the possible cause of the re¬ 
duced perpendicular diffusion in helical rotating flows. 
As the perpendicular energy spectrum in this case is 
steeper than in the absence of helicity, then the smaller 
energy at small scales results in less mixing and diffusion. 


2. Vertical diffusion 


Figure \T^ shows the mean vertical passive scalar con¬ 
centration 0{z), the mean vertical flux 6vz{z), and the 
pointwise value of Bziz) at different times in run 682 . In 
this case, the profiles are more similar to those obtained 
in the isotropic and homogeneous case: 6{z) and 0vz{z) 
are respectively symmetric and antisymmetric with re¬ 
spect to Z = TT. 

As in the case of horizontal diffusion, we can obtain 
the vertical turbulent diffusion coefficient as a function 
of time by computing the mean value of T>z{z,t) for all 
values of z. Figure 15 shows 'Dz{t) for runs ASj and BSz 
(both with Ro = 0.01, respectively without and with he¬ 
licity). Note that horizontal turbulent diffusion is larger 
in the presence of helicity, even more than in the isotropic 
case. 


D. Spatial distribution and structures 

Results shown above suggest that both horizontal and 
vertical diffusions are affected by rotation and by the 


FIG. 17: Passive scalar concentration in a vertical slice of 
runs A3z (left) and B3z (right) at time t = 0.5 (top row), and 
at time t — 1.5 (bottom row). 


presence of helicity. Figure shows a horizontal slice of 
the passive scalar concentration in runs A3x and BSa, at 
t = 1 (i.e., around the time the turbulent diffusion coef¬ 
ficients TAx and TAz reach a turbulent steady value). As 
also observed in [H], the initial Gaussian profile in the 
non-helical rotating flow (run AS^,) diffuses in time, and 
also bends and rotates. As previously mentioned, the 
overturning of the profile is caused by the Coriolis force 
(see also M)- In the helical rotating flow (run BS^,), we 
also observe this overturning, although the initial profile 
is less diffused (as indicated, e.g., by the most extreme 
values in the x-coordinate for which a significant concen¬ 
tration of the passive scalar can be observed, which are 
larger in run AS^,). 


Diffusion in the parallel direction in rotating flows is 


of a different nature than vertical diffusion (see Fig. 17 


which shows vertical slices of the passive scalar concen¬ 
tration in runs A3z and BSz at t = 0.5 and t = 1.5.). 
In the rotating non-helical case, the passive scalar initial 
profile is diffused in vertical stripes, created by updrafts 
or downdrafts inside columnar structures of the veloc¬ 
ity field El- These columnar structures in the velocity 
and vorticity fields have been reported in rapidly rotating 
flows, and are associated with the bidimensionalization of 
the flow [151t48) . As time increases, the stripes observed 
in the passive scalar in Fig. are further streched, re¬ 
sulting in larger mixing and diffusion. Note however that 
in the presence of helicity, the stripes are increased even 
further, in good agrement with the increased diffusion 
in helical flows reported above. This can be understood 
as the presence of helicity in the flow requires the three 
components of the velocity to be non-zero, resulting in a 
more three-dimensional flow. 
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V. CONCLUSIONS 

We analyzed data from direct numerical simulations 
of advection and diffusion of a passive scalar in rotating 
helical and non-helical turbulent flows. A total of 18 sim¬ 
ulations with spatial resolution of 512^ grid points was 
performed, using different Reynolds and Rossby num¬ 
bers, and changing the forcing and initial conditions of 
the passive scalar, to meassure energy and passive scalar 
spectra, anisotropic velocity and passive scalar structure 
functions, probability density functions, and diffusion co¬ 
efficients in the directions parallel and perpendicular to 
the rotations axis. 

In the first part of the paper we studied scaling laws of 
the energy and passive scalar variance, using spectra and 
structure functions in the horizontal and vertical direc¬ 
tions. We showed that helicity affects the inertial range 
scaling of the passive scalar, with its variance following 
a spectral law consistent with ^ This scaling is 

shallower than the one found for passive scalars in non¬ 
helical rotating turbulence [53], and consistent with a 
phenomenological argument that states that if the en¬ 
ergy follows a power law ^ in the inertial range, 
then the passive scalar variance should follow a power 
law ^ fc”"® with ne = (5 —n)/2. This argument, already 
proposed in [23j for rotating and non rotating non-helical 
flows, was found here to uphold also in the presence of 
helicity. The study of structure functions confirms these 
scaling laws, and indicates that the passive scalar is more 
anisotropic at small scales than velocity field. Also, the 
passive scalar was found to be more intermittent than 
the velocity field, a well known result, what which be¬ 


comes more pronounced in the presence of rotation and 
of helicity. The anomalous scaling exponents for the pas¬ 
sive scalar can be approximated using Kraichnan’s model 
with the second order exponent ^2 obtained from our phe¬ 
nomenological model, and with a dimensionality d = 2. 
As in the case of rotating non-helical flows studied previ¬ 
ously in [T], this value of d was interpreted as a result of 
the quasi-bidimentionalization of the distribution of the 
passive scalar in the presence of rotation. 

In the second part of the paper, the analisys of the 
effective diffusion coeficients calculated from Pick’s law 
show that for isotropic flows (i.e., without rotation) helic¬ 
ity increases turbulent diffusion, in agreement with pre¬ 
vious models and theoretical predictions [iiiiinj- In the 
presence of rotation, results indicate that the overall ef¬ 
fect of rotation (irrespectively of the content of helicity 
of the flow) is to decrease horizontal diffusion, while the 
effect on vertical diffusion is less pronounced. Helicity 
further decreases horizontal diffusion but increases ver¬ 
tical diffusion (compared with the non-helical rotating 
case). The decrease in horizontal diffusion was explained 
with a simple model for turbulence diffusivity based on 
the available energy for the small-scale turbulent fluctu¬ 
ations. 
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